Biological Evaluation and Molecular Docking Studies of Novel 1,3,4-Oxadiazole Derivatives of 4,6-Dimethyl-2-sulfanylpyridine-3-carboxamide

To date, chronic inflammation is involved in most main human pathologies such as cancer, and autoimmune, cardiovascular or neurodegenerative disorders. Studies suggest that different prostanoids, especially prostaglandin E2, and their own synthase (cyclooxygenase enzyme-COX) can promote tumor growth by activating signaling pathways which control cell proliferation, migration, apoptosis, and angiogenesis. Non-steroidal anti-inflammatory drugs (NSAIDs) are used, alongside corticosteroids, to treat inflammatory symptoms particularly in all chronic diseases. However, their toxicity from COX inhibition and the suppression of physiologically important prostaglandins limits their use. Therefore, in continuation of our efforts in the development of potent, safe, non-toxic chemopreventive compounds, we report herein the design, synthesis, biological evaluation of new series of Schiff base-type hybrid compounds containing differently substituted N-acyl hydrazone moieties, 1,3,4-oxadiazole ring, and 4,6-dimethylpyridine core. The anti-COX-1/COX-2, antioxidant and anticancer activities were studied. Schiff base 13, containing 2-bromobenzylidene residue inhibited the activity of both isoenzymes, COX-1 and COX-2 at a lower concentration than standard drugs, and its COX-2/COX-1 selectivity ratio was similar to meloxicam. Furthermore, the results of cytotoxicity assay indicated that all of the tested compounds exhibited potent anti-cancer activity against A549, MCF-7, LoVo, and LoVo/Dx cell lines, compared with piroxicam and meloxicam. Moreover, our experimental study was supported by density functional theory (DFT) and molecular docking to describe the binding mode of new structures to cyclooxygenase.

In our previous paper, we have reported the synthesis of twelve different Schiff base derivatives of N-(2-hydrazinyl-2-oxoethyl)-4,6-dimethyl-2-sulfanylpyridine-3carboxamide ( Figure 1) [55]. The compounds were tested for anti-COX, antioxidant and anticancer activity. The derivatives showed promising properties, therefore we decided to obtain analogous structures additionally containing the pharmacophoric ring of 1,3,4oxadiazole. The heterocyclic moiety may enhance the pharmacological action of the entire molecule.
Based on the above information, it can be concluded that our newly synthesized derivatives, which are hybrid compounds connecting the structure of 4,6dimethylpyridine, 1,3,4-oxadiazole ring and N-acyl hydrazone moiety, may exhibit significant anti-inflammatory and chemopreventive effects, and their mechanism of action may be associated with anti-COX, antioxidant and anticancer activities ( Figure 1). Considering that the binding pocket of COX-2 isoenzyme is bigger than that of COX-1, we expect that the extensive structure of new molecules will enhance the COX-2 selectivity compared to the previous Schiff bases, which inhibited COX-1 at lower concentrations than COX-2 [56]. Based on the above information, it can be concluded that our newly synthesized derivatives, which are hybrid compounds connecting the structure of 4,6-dimethylpyridine, 1,3,4oxadiazole ring and N-acyl hydrazone moiety, may exhibit significant anti-inflammatory and chemopreventive effects, and their mechanism of action may be associated with anti-COX, antioxidant and anticancer activities ( Figure 1). Considering that the binding pocket of COX-2 isoenzyme is bigger than that of COX-1, we expect that the extensive structure of new molecules will enhance the COX-2 selectivity compared to the previous Schiff bases, which inhibited COX-1 at lower concentrations than COX-2 [56].
The new compounds were examined as prospective COX-1/COX-2 inhibitors. Moreover, antioxidant and cytotoxic properties against A549 (pulmonary basial cell alveolar adenocarcinoma), MCF-7 (breast adenocarcinoma), LoVo (colon adenocarcinoma), and its drug resistant subline LoVo/Dx cell lines were studied to check the chemopreventive potential of the compounds. In addition, the structure study was performed on the basis of density functional theory (DFT) to investigate the energetic properties of the studied compounds depending on their conformation. Finally, molecular docking was carried out in order to suggest the binding mode and correlate it with biological activity.

Chemistry
The synthesis of N-(2-hydrazinyl-2-oxoethyl)-4,6-dimethyl-2-sulfanylpyridine-3-carboxamide 1 was performed according to the protocols published previously [19]. Scheme 1 presents the synthesis of compounds which have not been described in the literature yet. Analytical and spectroscopic properties of all newly obtained derivatives were in good agreement with their predicted structures and are summarized in the experimental section. Initially, compound 1 was heated at reflux in the presence of carbon disulfide in basic conditions in ethanol. Subsequently, the reaction mixture was poured onto crushed ice and acidified with hydrochloric acid. As a result, hydrazide 1 was subjected to intramolecular cyclization and formation of a five-membered 1,3,4-oxadiazol-2-thione ring. In the next step, compound 2 was reacted with ethyl bromoacetate in a water-ethanolic potassium hydroxide solution. Thus, the ester derivative 3 was obtained. Subsequently, compound 3 was converted into hydrazide 4 by reaction with hydrazine hydrate. The synthesis of nine different Schiff base derivatives of compound 4 was carried out by the reaction of compound 4 with several aromatic aldehydes in methanol in the presence of a catalytic amount of acetic acid with good yields. The purity of the synthesized compounds was checked by elemental analyses. The structures of the various synthesized compounds were determined based on spectral data analysis, such as FT-IR and 1 H NMR and 13 C NMR. The new compounds were examined as prospective COX-1/COX-2 inhibitors. Moreover, antioxidant and cytotoxic properties against A549 (pulmonary basial cell alveolar adenocarcinoma), MCF-7 (breast adenocarcinoma), LoVo (colon adenocarcinoma), and its drug resistant subline LoVo/Dx cell lines were studied to check the chemopreventive potential of the compounds. In addition, the structure study was performed on the basis of density functional theory (DFT) to investigate the energetic properties of the studied compounds depending on their conformation. Finally, molecular docking was carried out in order to suggest the binding mode and correlate it with biological activity.

Chemistry
The synthesis of N-(2-hydrazinyl-2-oxoethyl)-4,6-dimethyl-2-sulfanylpyridine-3-carboxamide 1 was performed according to the protocols published previously [19]. Scheme 1 presents the synthesis of compounds which have not been described in the literature yet. Analytical and spectroscopic properties of all newly obtained derivatives were in good agreement with their predicted structures and are summarized in the experimental section. Initially, compound 1 was heated at reflux in the presence of carbon disulfide in basic conditions in ethanol. Subsequently, the reaction mixture was poured onto crushed ice and acidified with hydrochloric acid. As a result, hydrazide 1 was subjected to intramolecular cyclization and formation of a five-membered 1,3,4-oxadiazol-2-thione ring. In the next step, compound 2 was reacted with ethyl bromoacetate in a water-ethanolic potassium hydroxide solution. Thus, the ester derivative 3 was obtained. Subsequently, compound 3 was converted into hydrazide 4 by reaction with hydrazine hydrate. The synthesis of nine different Schiff base derivatives of compound 4 was carried out by the reaction of compound 4 with several aromatic aldehydes in methanol in the presence of a catalytic amount of acetic acid with good yields. The purity of the synthesized compounds was checked by elemental analyses. The structures of the various synthesized compounds were determined based on spectral data analysis, such as FT-IR and 1 H NMR and 13 C NMR. The FT-IR spectrum of hydrazide 4 showed peaks at 1648 and 1623 cm −1 due to two carbonyl functions derived from the amide and hydrazide structure. Additionally, the IR spectra of compounds 4 and 5-13 exhibited in the 3294-3035 cm −1 range, the NH 2 and NH weak band of the NHNH 2 , CONH and CONH-N= functions.
The 1 H NMR spectrum of hydrazide 4 displayed no signals belonging to the OCH 2 CH 3 group; instead, new signals derived from the hydrazide structure appeared at 6.07 (NHNH 2 ) and 9.29 (NHNH 2 ) ppm integrating for two protons and one proton, respectively. The 1 H NMR spectra of compounds 5-13 displayed additional signals due to the azomethine group and the aromatic ring derived from the aldehyde moiety at the aromatic region, while the signal belonging to NH 2 group of the hydrazide structure did not appear.
It should be noted that the 1 H NMR spectra of compounds 5-13 show double signals corresponding to two amidic NH protons, azometine moiety protons and methylene protons. According to the literature, compounds having an arylidene-hydrazide structure may exist as E/Z geometrical isomers around the C=N double bond and as cis/trans amide conformers. It has been reported that when compounds containing an imine bond are dissolved in dimethyl-d 6 sulfoxide, they are present in the form of a geometrical E isomer. The Z isomer can be stabilized in less polar solvents by an intramolecular hydrogen bond [57][58][59][60]. In this study, the spectral data were obtained in dimethyl-d 6 sulfoxide solution and no signal belonging to Z isomer was observed. On the other hand, the cis-trans conformers of E isomer were present in the dimethyl-d 6 sulfoxide solution of compounds 5-13. In the 1 H NMR spectra of 5-13, in particular, two sets of signals each belonging to the CH 2 group, N=CH group and CONH group of cis and trans conformers were observed between 3.96 and 4.44, 7.96 and 8.54, 11.56 and 11.99 ppm respectively. The upfield lines of protons of the listed groups were assigned to the cis-conformer of the amide structure, while the downfield lines of the protons of the same groups were assigned to the trans-conformer of the amide structure [59].

Cyclooxygenase Inhibition Assay
The compounds were studied for their potencies to inhibit COX-1 and COX-2 enzymes by the colorimetric inhibitor screening assay. The IC 50 values (i.e., the concentration of tested compounds (µM) that can exert 50% inhibition of the enzyme activity) and the COX-2/COX-1 selectivity ratios were calculated after 2 min of incubation for each investigated and reference compound (piroxicam and meloxicam).
As can be seen from the data presented in Table 1, hydrazide 4 and its 4methylsulfanylbenzylidene derivative 8 showed only COX-1 inhibitory activity, while the compound with 4-chlorobenzylidene moiety 10 selectively inhibited the activity of COX-2 isoenzyme. Four of the Schiff bases 5, 11, 12 and 13 with unsubstituted benzylidene moiety or containing 4-trifluoro, 3-chloro, 2-bromo groups, respectively exhibited higher anti-COX-1 activity compared with piroxicam and meloxicam. What is more, compound 13 inhibited the activity of COX-2 at a lower concentration than standards, and its COX-2/COX-1 selectivity ratio was similar to meloxicam.

MTT Assay
The MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl tetrazolium bromide) assay was used to determine the cytotoxicity of tested compounds towards different cell lines: normal human dermal fibroblasts (NHDF), kidney epithelial cells (VERO), fibroblasts from Chinese hamster lung (V79) and four human cancer cells: A549 (pulmonary basal cell alveolar adenocarcinoma), MCF-7 (breast adenocarcinoma), LoVo (colon adenocarcinoma), and its drug resistant subline LoVo/Dx. The cell lines used in this study express COX-1 and COX-2. The LoVo line expressed COX-1, the LoVo/Dx subline express both COX-1 and COX-2 [61]. The A549 and MCF-7 lines shows expression of both enzymes [62,63]. This assay is suitable for the measurement of drug sensitivity in cultured cells, specified as the concentration of the compound required to achieve 50% growth inhibition, compared to the growth of the control cultured without any drug (50% inhibitory concentration, IC 50 ) [64]. The results of the experiment were made for screening purpose, and are presented in Table 2. Note: Data are shown as a mean ± standard deviation, NA stands for "not applicable".
Firstly, the effect of the new compounds on healthy cells (NHDF, V79 and VERO) were performed. With regards to NHDF cells, both, hydrazide 4 and all of the Schiff bases exhibited similar toxicity (IC 50 in the range of 57.47-59.9 µM). The VERO line exhibited more differences in cells sensitivity: there was no toxicity caused to VERO cells by standard compounds (piroxicam and meloxicam) and IC 50 for tested compounds was in the range 50.38 µM for compound 8 to 129.8 µM for compound 10. The results for V79 cell line were more similar to NHDF cells: IC 50 was in the range 47.70 µM (for compound 7) to 84.30 µM (for compound 5).
The anti-cancer activity of the tested compounds was established towards A549 cell line. A549 cells are adenocarcinomic human alveolar basal epithelial cells, used as a model of lung adenocarcinoma [65]. Inhibitory concentrations caused 50% of growth inhibition for these cells were similar for all tested compounds (3.16-4.34 µM), and was 25 to 35 times lower than that of piroxicam and meloxicam. The next tested cell line was MCF-7 from human breast adenocarcinoma. In this case, the standard compounds possessed no antitumor activity in the range of used concentrations, while the inhibitory concentrations of tested compounds were in the range of 6.33-9.63 µM. Antitumor activity towards LoVo cells (from human colon) was similar as for A549 cells, and the inhibitory concentrations were 34-to 45-fold lower than that of piroxicam and meloxicam. The drug resistant subline LoVo/Dx exhibited more different sensitivity to the tested compounds. Hydrazide 4 and five Schiff bases 5, 6, 7, 9 and 13 with unsubstituted benzylidene moiety or containing 4-fluoro, 4-methyl, 4-cyano and 2-bromo groups, respectively inhibited the growth of LoVo/Dx cells in the concentrations 13-to 22-fold lower than the standards.
The performed investigation revealed that all tested compounds had more significant anticancer activity than piroxicam and meloxicam. Moreover, the therapeutic index, i.e., the ratio of the concentrations that inhibit 50% of healthy and cancerous cells (A549, MCF-7 and LoVo cell lines) was high for all investigated compound, and the IC 50 values for tumor cells were 1.3-and even 27-fold lower than for normal cells (Table 3). Note: Data are shown as a mean ± standard deviation. a Represents the significance level at p < 0.05.

Evaluation of Reactive Oxygen Species (ROS) Level Inside the Cells
ROS participates in various redox-regulatory mechanisms of cells to maintain homeostasis. When the system's ability to neutralize and eliminate free radicals and active intermediates fails, the intracellular redox potential shifts towards oxidative stress leading to DNA mutations, which can further promote neoplastic transformation. A chronic inflammation is the main system capable of inducing an oxidative stress. Prostaglandins induce the expression of certain inflammatory cytokines, which can, in turn, enhance the production of ROS. The inhibition of these paths at an early stage of neoplasia may be a critical moment in the process of chemoprevention. [66]. Taking this into account, the next part of the investigations was the ROS scavenging activity of studied compounds examined under or without the influence of oxidative stress, generated by the presence of H 2 O 2 inside the cell cultures. As in the previous tests, two oxicams (meloxicam and piroxicam) were taken as the reference drugs (Tables 4 and 5) due to their proven ability to scavenging ROS [67][68][69].
Two cell lines were chosen for investigations: V79 (Chinese hamster lung fibroblasts) and LoVo. V79 normal cells are widely used as the cell line tested in the studies of the oxidative stress [70][71][72]. These cells are a living system model to explain the mechanism of ROS effects caused by varied compounds [73].
The tests established on V79 cell line revealed that among ten tested compounds, two of them (hydrazide 4 and Schiff base 13 containing 2-bromobenzylidene group) showed significant ROS scavenging activity under normal conditions almost two times higher than the standard drugs. The same compounds exhibited the best values during oxidative stress, induced by H 2 O 2 ( Table 4).  Considering the test on LoVo line 4-methylsulphanylbenzylidene derivative (8) showed high ROS scavenging activity under normal conditions. On the other hand the best activity during oxidative stress characterized compounds 4 and 13 as in the case of the first cell line V79 (Table 5, Figure 2).
Comparing the antioxidant activity of our ten tested compounds to standard ROS scavengers: trolox and ascorbic acid, we can see that the activity of all compounds tested on V79 line, is lower than trolox and ascorbic acid (Table 4). Trolox and ascorbic acid antioxidant activities were previously tested by Szczęśniak-Sięga et al. on V79 cell line in the similar conditions, and the values of their mean E/E 0 were similar to our results [74]. The same situation is observed for LoVo line, where the antioxidant activity of tested compounds were also lower than for standard scavengers, but very similar for all tested compounds ( Table 5).
The same situation is observed for LoVo line, where the antioxidant activity of tested compounds were also lower than for standard scavengers, but very similar for all tested compounds (Table 5).
According to the results, it can be seen that our tested compounds can protect healthy cells against oxidative stress to a greater extent than cancer cells. Results are presented as RFU (relative fluorescent units) of fluorescence levels. The statistical significance of the differences between the results for the tested compounds, compared to the control, was calculated using Tukey's post hoc test (* < 0.05).

Static Density Functional Theory (DFT) Models and Molecular Docking Study
Molecular modeling techniques are useful in development of the interaction models between ligands and receptors. Within this study we applied diverse methods to characterize the receptors (sequence similarity evaluation), the ligands themselves (density functional theory, DFT) and their interactions with protein receptors (docking studies). We will start with the discussion of the protein models. For human COX-1 (PDB-Protein Data Bank code: 6Y3C [75]), the sequence similarity to ovine protein model (PDB code: 4O1Z According to the results, it can be seen that our tested compounds can protect healthy cells against oxidative stress to a greater extent than cancer cells.

Static Density Functional Theory (DFT) Models and Molecular Docking Study
Molecular modeling techniques are useful in development of the interaction models between ligands and receptors. Within this study we applied diverse methods to characterize the receptors (sequence similarity evaluation), the ligands themselves (density functional theory, DFT) and their interactions with protein receptors (docking studies). We will start with the discussion of the protein models. For human COX-1 (PDB-Protein Data Bank code: 6Y3C [75]), the sequence similarity to ovine protein model (PDB code: 4O1Z [76]) was equal to 91% and for the human COX-2 (PDB code: 5KIR [77]) the resemblance to the murine model (PDB code: 4M11 [76]) reached 88%-the sequence alignments are presented in Figure S1 of Supplementary Material. Sequence homology between human receptors was close to 66%, hence the different mode of binding, as well as affinity energies, were anticipated. Significant homology in the binding places between COX-1 orthologes was seen-the binding site of the protein structure with 6Y3C code (apo form, with no ligand molecule) was inferred from knowledge of the 4O1Z binding pocket-molecular docking has shown that our assumption was correct and even if affinity energies for human COX-1 were slightly lower, the docking procedure itself was successful. Residues such as Leu352, PhE-518, Ser530, Trp387 and Val349 were involved in the direct interactions with ligands in all the studied receptors.
Two types of conformations of the studied set of compounds were found in the DFT structural optimizations: one that assumed more packed, sandwich-like structure-denoted by additional "S" in brackets after the compound number-and the second one which was more extended (with "E" in brackets after the compound number). The details are shown in the Figure 3 and in the Supplementary Materials (Table S1). to the murine model (PDB code: 4M11 [76]) reached 88%-the sequence alignments are presented in Figure S1 of Supplementary Material. Sequence homology between human receptors was close to 66%, hence the different mode of binding, as well as affinity energies, were anticipated. Significant homology in the binding places between COX-1 orthologes was seen-the binding site of the protein structure with 6Y3C code (apo form, with no ligand molecule) was inferred from knowledge of the 4O1Z binding pocket-molecular docking has shown that our assumption was correct and even if affinity energies for human COX-1 were slightly lower, the docking procedure itself was successful. Residues such as Leu352, PhE-518, Ser530, Trp387 and Val349 were involved in the direct interactions with ligands in all the studied receptors.
Two types of conformations of the studied set of compounds were found in the DFT structural optimizations: one that assumed more packed, sandwich-like structure-denoted by additional "S" in brackets after the compound number-and the second one which was more extended (with "E" in brackets after the compound number). The details are shown in the Figure 3 and in the Supplementary Materials (Table S1). The structures exhibited diverse features, including intramolecular N-H...O hydrogen bonds and non-covalent stacking-like interactions between the phenyl and pyridine rings, which affected the located conformations and their stability. The set of compounds shares the common part and differs generally only by the phenyl ring substitution. The applied energetic criterion-relative energy between the conformers of a given compound-shows that the sandwich-like structures are in general ca. 8-10 kcal/mol lower in energy. This energy separation between conformers was found in the gas phase DFT The structures exhibited diverse features, including intramolecular N-H . . . O hydrogen bonds and non-covalent stacking-like interactions between the phenyl and pyridine rings, which affected the located conformations and their stability. The set of compounds shares the common part and differs generally only by the phenyl ring substitution. The applied energetic criterion-relative energy between the conformers of a given compoundshows that the sandwich-like structures are in general ca. 8-10 kcal/mol lower in energy. This energy separation between conformers was found in the gas phase DFT study based on different initial arrangement of the central 1,3,4-oxadiazole ring. Both types of conformations were then examined via the docking approach.
In the next step, the protein-ligand complex was analyzed using the flexible docking method. The idea of the flexible docking is that both the ligand and the selected receptor residues are treated as conformationally flexible. The initial conformation of the ligand is modified to adapt to the binding site and the selected residues of the binding site are also allowed to reorient. Initial docking validation was performed, and it consisted of re-docking of the ligands that were previously co-crystallized with the corresponding receptors. It was meloxicam (MXM) for both ovine COX-1 (4O1Z) and murine COX-2 (4M11), and rofecoxib (RCX) for COX-2 (5KIR) from Homo sapiens. For human COX-1 (6Y3C) validation was not possible because of lack of the co-crystallized ligand alongside the protein in the PDB repository. The validation set of the docked structures was very similar to the crystallized cases (see Figure 4): the RMSD (root-mean-square deviation) factor for the ligand position was equal to 1.418 Å for rofecoxib in 5KIR, 0.930 Å for meloxicam in 4M11 and 1.037 Å for meloxicam in 4O1Z. The RMSD was calculated based on the structures with the highest conformational similarity to the co-crystallized ligands. For every structure of the MXM and RCX, the lowest RMSD value was obtained for conformations with the second lowest affinity energy.
method. The idea of the flexible docking is that both the ligand and the selected receptor residues are treated as conformationally flexible. The initial conformation of the ligand is modified to adapt to the binding site and the selected residues of the binding site are also allowed to reorient. Initial docking validation was performed, and it consisted of re-docking of the ligands that were previously co-crystallized with the corresponding receptors. It was meloxicam (MXM) for both ovine COX-1 (4O1Z) and murine COX-2 (4M11), and rofecoxib (RCX) for COX-2 (5KIR) from Homo sapiens. For human COX-1 (6Y3C) validation was not possible because of lack of the co-crystallized ligand alongside the protein in the PDB repository. The validation set of the docked structures was very similar to the crystallized cases (see Figure 4): the RMSD (root-mean-square deviation) factor for the ligand position was equal to 1.418 Å for rofecoxib in 5KIR, 0.930 Å for meloxicam in 4M11 and 1.037 Å for meloxicam in 4O1Z. The RMSD was calculated based on the structures with the highest conformational similarity to the co-crystallized ligands. For every structure of the MXM and RCX, the lowest RMSD value was obtained for conformations with the second lowest affinity energy. After successful validation of the docking protocol, the further characterization of the binding pocket interactions with the set of the studied molecules was performed. The complicated network of interactions, difficult to present in a 3D representation is provided in the 2D form in Figure S2 of the Supplementary Material. The most favorable structures with regards to the binding affinity energy (presented in Figures 5 and 6) were chosen for After successful validation of the docking protocol, the further characterization of the binding pocket interactions with the set of the studied molecules was performed. The complicated network of interactions, difficult to present in a 3D representation is provided in the 2D form in Figure S2 of the Supplementary Material. The most favorable structures with regards to the binding affinity energy (presented in Figures 5 and 6) were chosen for detailed analysis. The most active extended-conformation compound towards the human COX-1 was 11(E). The stability of the binding of this compound by the human COX-1 can be explained by formation of hydrogen bonds between the Arg120, Ser530 and Asn375 residues and the ligand (depicted in Figure S2). Stabilizing non-covalent interaction (halogen bonding) was detected by the Discovery Studio Visualizer 2021 [78] between one of the fluorine substituents of the phenyl ring and the Phe-529 residue (also see Figure S2). Other non-covalent interactions present in the binding pocket are graphically presented in Figure S2. The structure represented by 11(S) compound was also the most stabilized among the others of sandwich-like conformations. In this case, the most important modes of the binding contained also the hydrogen bonds between the ligand and the Arg120, Asn375 residues. Numerous other non-covalent interactions were detected, and they are presented in Figure S2.
be explained by formation of hydrogen bonds between the Arg120, Ser530 and Asn375 residues and the ligand (depicted in Figure S2). Stabilizing non-covalent interaction (halogen bonding) was detected by the Discovery Studio Visualizer 2021 [78] between one of the fluorine substituents of the phenyl ring and the Phe-529 residue (also see Figure S2). Other non-covalent interactions present in the binding pocket are graphically presented in Figure S2. The structure represented by 11(S) compound was also the most stabilized among the others of sandwich-like conformations. In this case, the most important modes of the binding contained also the hydrogen bonds between the ligand and the Arg120, Asn375 residues. Numerous other non-covalent interactions were detected, and they are presented in Figure S2.  It is worth noting that the docking to the apo structure of the human COX-1 was successful because of the use of the flexible docking protocol and previous characterization of the similarity of the binding sites with the use of bioinformatics tools (see Figure S1). For the ovine COX-1, the affinity energies were generally larger by ca. 2 kcal/mol for each compound comparing with the previous receptor (Tables 6 and 7). The most selective extended structure was 7(E) and the best among the sandwich-like structures was 9(S). For the 7(E) compound a large number of π-π stacked, amide-π stacked, π-sulphur and π-σ interactions were present (see Figure S2). In the case of 9(S) the π-alkyl, π-σ, alkyl and also π-sulphur interactions were the most significant. The binding energy for both 7(E) and 9(S) compounds exceeded -10.0 kcal/mol. The extended 11(E) compound was found as the most active for the Mus musculus COX-2 ( Table 8). The halogen bonds formation between fluorine substituents of ligand phenyl ring and Met522 residue was observed. The Arg120 residue and the nitrogen from the ligand pyridine ring were involved in the hydrogen bonding; non-covalent interactions are depicted in Figure S2. Its affinity towards Mus musculus COX-2 was even higher than that of the co-crystallized ligand-MXM. Among sandwich-like structures the most promising seems to be the 9(S) molecule. In this case, the hydrogen bonds between the ligand and the residues Ser530 and Arg120 are present. Plethora of Van der Waals interactions between alkyl and aromatic parts of ligand and valine and leucine residues were also present (see Figure S2). The human COX-2 was bounded more efficiently by extended conformation of the examined compounds ( Table 9). The differences were most significant for 11(E) where the affinity energy between two investigated structures varied by more than 2 kcal/mol. The 11(E) molecule in the binding pocket was stabilized by the network of intermolecular hydrogen bonds between the ligand and the Thr94, His90, Ser353, Tyr355 and Arg120 residues. Numerous non-covalent interactions were detected, and they are graphically presented in Figure S2. The 13(S) molecule had the highest affinity towards the receptor among the sandwich-like structures. The major role was played by the amide-π stacked and π-σ as well as other non-covalent interactions present in the binding pocket (for details see Figure S2). It is worth noting that the docking to the apo structure of the human COX-1 was successful because of the use of the flexible docking protocol and previous characterization of the similarity of the binding sites with the use of bioinformatics tools (see Figure  S1). For the ovine COX-1, the affinity energies were generally larger by ca. 2 kcal/mol for each compound comparing with the previous receptor (Tables 6 and 7). The most selective extended structure was 7(E) and the best among the sandwich-like structures was 9(S).     Concluding the docking study, we have found the substituent effect to change the affinity energy by 2-3 kcal/mol, which is a moderate but visible effect (ca. 20%). It was also found that the sandwich-like structures varied less in the affinity energies than the extended conformations. This is especially for the human COX-1, where only molecule conformations for structure 4 possessed lower affinity comparing to other compounds. The affinity energy differences between the most and the least favorable conformations of each of the ligands were ca. 3-4 kcal/mol. The affinities for the most stable structures varied even less. Generally speaking, the lower affinity energies were obtained when the docking procedure was applied to human and murine COX-2 receptors. The compound 4 was the least selective for the examined receptors, probably due to the lack of stabilizing non-covalent interactions provided by the phenyl ring and its substituents. It was most pronounced for receptors with 6Y3C, 4M11 and 5KIR PDB codes. The 11(S) compound could be considered as potentially the most promising structure and it seems that the CF 3 substituent of the phenyl ring plays an important role in stabilizing that ligand in the binding pocket (due to the intermolecular halogen bonds formation).

Instruments and Chemicals
All solvents, reagents and chemicals used during experiments described in this paper were delivered by commercial suppliers (Alchem, Wrocław, Poland; Chemat, Gdańsk, Poland; Archem, Łany, Poland) and were used without further purification. Any dry solvents were received due to standard procedures. Reaction progress was monitored by the thin-layer chromatography (TLC) technique, on TLC plates made of 60-254 silica gel, and was visualised by UV light at 254/366 nm. Melting points of final compounds were determined on Electrothermal Mel-Temp 1101D apparatus (Cole-Parmer, Vernon Hills, IL, USA) using open capillary method, no correction needed. 1 H NMR (300 MHz) and 13 C NMR (75 MHz) spectra were recorded using Bruker 300 MHz NMR spectrometer (Bruker Analytische Messtechnik GmbH, Rheinstetten, Germany) in DMSO-d 6 , with tetramethylsilane (TMS) as an internal reference. Chemical shifts (δ) were reported in ppm. In order to record and read spectra, TopSpin 3.6.2. (Bruker Daltonik, GmbH, Bremen, Germany) program was used. Elemental analyses for carbon, nitrogen and hydrogen were carried out on a Carlo Erba NA 1500 analyzer and were within ±0.4% of the theoretical value. FT-IR spectra were measured on Nicolet iS50 FT-IR Spectrometer (Thermo Fisher Scientific, Waltham, MA, USA). Frequencies were reported in cm −1 . All samples were solid, and spectra were read by OMNIC Spectra 2.0 (Thermo Fisher Scientific, Waltham, MA, USA).

Preparation and Experimental Properties of Compounds 2-13
The synthesis protocols and experimental data for compound 1 have already been reported [19].
The hydrazide 1 (1.27 g, 0.005 mol) and KOH (0.56 g, 0.01 mol) were dissolved in ethanol (50 mL) in a round bottom flask. To this stirring mixture carbon disulphide (3 mL, 0.05 mol) was added and the whole was refluxed for 10 h till evolution of hydrogen sulfide was ceased. Then, the reaction mixture was cooled and slowly acidified with diluted hydrochloric acid. Formed precipitate was filtered off, washed with cold water, dried and recrystallized from ethanol giving white solid of compound 2.
Yield In a 100 mL round bottom flask, 3.68 g (0.01 mol) of compound 3 and 30 mL of methanol were placed. The resulting mixture was refluxed with stirring. After the compound 3 had dissolved, 5.1 mL of hydrazine hydrate was added, and the mixture was heated for 5 h. The obtained precipitate was filtered off and allowed to dry. The compound 4 was then recrystallized in methanol. Yield

Cell Lines
Six cell lines were used in the investigations. Three of them were normal: NHDF (normal human dermal fibroblasts), purchased from Lonza (Verviers, Belgium), V79 (fibroblasts from Chinese hamster lung) and VERO (kidney epithelial cells) obtained from ECACC (European Collection of Authenticated Cell Cultures). The rest of cells were human cancer cell lines: LoVo (colon adenocarcinoma), their drug resistant subline LoVo/Dx and A549 (pulmonary basial cell alveolar adenocarcinoma) cell line. Cancer lines were also obtained from ECACC.

Cell Culture Conditions
Each cell line was grown in the culture media recommended by supplier. Before the test they were detached with trypsin/EDTA solution, then to neutralize the effect of Trypsin/EDTA solution, FBS containing medium was used. After centrifugation the cells were stained with 0.4% solution of trypan blue, counted and inspected for viability using microscope. In the end, cells were inserted into 96-well culture plates and incubated in CO 2 incubator (37 • C, 24 h). The number of cells was 2 × 10 3 cells per well. The tested compounds with different concentrations (5, 10, 20, 50, 100 µM) were dissolved in DMSO and then added to the cells (the final DMSO concentration was 0.1%). The cultures were incubated for another 48 h. After that time the cells were collected to be used in tests.

Cyclooxygenase Inhibitory Activity
A COX Colorimetric Inhibitor Screening Assay Kit, produced by Cayman Chemical Company, Ann Arbor, MI, USA, was used. It used the colorimetric monitoring of oxidized form of TMPD (N,N,N ,N -tetramethyl-p-phenylenediamine) which was produced during reduction of prostaglandin G 2 (PGG 2 ) to PGH 2 . The change of colour was observed and measured spectrofotometrically at 590 nm. Reagents used it this assay were: COX-1 and COX-2 enzymes, Tris-HCl buffer, solutions of heme in DMSO and TMPD, arachidonic acid, KOH. To be sure that 100% enzymatic activity was achieved each sample was measured 3 times. The probes were measured after two minutes of incubation with tested compounds in the comparison to the initial activity of enzyme. It allowed to determine IC 50 value, where 50% inhibition of the enzyme activity was observed.

MTT Assay
MTT assay was used to find out how tested compounds influence the metabolic activity of investigated cell lines. The cells were incubated with tested compounds. After removing of supernatant, 1 mg/mL of MTT solution in MEM was added to the plate (to each well) and then the plates were incubated at 37 • C for 2 h. The medium was removed again, and formazan crystals were dissolved in isopropanol. A Varioscan LUX microplate reader (Thermo Fisher Scientific, Waltham, MA, USA) was used to measure absorbance at 570 nm.

Estimation of Intracellular ROS Level
First, the tested compounds (in concentration of 100 µM) were added to the cell cultures and cells were incubated for 4 h. Then cells were washed and incubated with DCFH-DA (non-fluorescent probe) in dark at 37 • C for 2 h at CO 2 [79]. Non-fluorescent, non-polar DCFH-DA (2 7 -dichlorodihydrofluorescein diacetate) at concentration of 25 µM was used as a marker of oxidative stress to determine the intracellular ROS levels. Then the cells were washed with PBS two times and 100 µL of H 2 O 2 were added for 30 min, which is proper time to decompose all H 2 O 2 [80]. In this time, DCFH-DA which penetrated into the cells, was hydrolyzed by esterases to polar, non-fluorescent DCFH (2 7dichlorodihydrofluorescein) and then oxidized in the presence of ROS (reactive oxygen species) to fluorescent DCF (2 7 -dichlorofluorescein). DCF was measured using Varioscan LUX microplate reader (Thermo Fisher Scientific, Waltham, MA, USA) at λ ex = 485 nm and λ em = 535 nm. Then results were presented as E/E 0 where E is test sample value and E 0 -control value.

Statistics
All results in Tables are presented as mean ± SEM (standard error of the mean) relative to the control (E/E 0 ). E is the culture with the tested compound and E 0 is the negative control without the compound. The routine statistical methods: two-way analysis of variance ANOVA as well as Tukey's post-hoc test were used. It allowed us to determine the statistical significance of results, where p < 0.05 was set as significant. The statistical indicators were estimated using Statistica 13.3 software.

Molecular Modeling-Computational Methodology
Initial models of the investigated compounds (see Figure 1) were constructed in the Molden 6.6 program [81] with the aim of quantum-chemical structural optimization before the main docking runs. The structure of each compound was optimized in accordance with default procedures in the Gaussian 16 suite of programs [82]. Energy minimization was carried out with the ωB97XD functional [83] derived within the framework of density functional theory (DFT) [84,85], with the correlation-consistent Dunning basis set denoted as cc-pVDZ [86]. In order to confirm that the optimized geometries correspond to the energy minima on the potential energy surface (PES), harmonic frequencies calculations were also conducted (yielding no imaginary frequencies). Resulting energy values were used to estimate the conformational preferences. Further, docking runs were carried out for the studied series of compounds. The structures of the receptors were taken from the Protein Data Bank (PDB) [87]. The compounds were docked into COX-1 (PDB code 4O1Z) [76] from Ovis aries and COX-2 (PDB code 4M11) [76] from Mus musculus models, which were loaded with the ligand (meloxicam, MXM). Additionally, COX-1 (PDB code 6Y3C) [75] and COX-2 (PDB code 5KIR) [77], cyclooxygenases counterparts from Homo sapiens, were also taken into consideration as receptors-of these, only the human COX-2 was loaded with rofecoxib (RCX) ligand, while the human COX-1 was in the apo form. Sequence alignments and similarity were generated with the ClustalW software [88]. Water molecules and cofactors that were in proximity to the binding site of receptors were removed. The Gasteiger charges and polar hydrogen atoms were added to the examined receptors and ligands. The flexible parts of the receptors were inferred from the knowledge of the center of cocrystallized ligands and the grid was set to the dimensions of 20 × 20 × 20 Å. The flexible part of the receptor was defined as residues located within 3.5 Å for COX-1 and within 4.0 Å for COX-2 from the centers of the binding pockets. The docking calculations were performed with exhaustiveness of 32 and energy range equal to 10 kcal/mol. The validation of the procedure was carried out with usage of the MXM and RCX ligands from the PDB database and, solely for the purpose to check the validation of the used computational setup, the polar hydrogens were not added. Editing of the protein structure (water and ligand removal) was carried out with the assistance of the VMD 1.9.3 program [89]. The preparation of the ligand and receptor, molecular docking and the visualization of the results were performed with the AutoDockTools 1.5.7 [90], AutoDock Vina 1.1.2 [91] and the open-source PyMOL 2.3.0 [92] packages. The 2D diagrams were prepared in the Discovery Studio Visualizer 2019 [78] with the use of the most energetically favorable geometries of the examined compounds, and the script vina_split of the AutoDock Vina 1.1.2 program [91] was used to split merged conformations from one PDBQT file into separate entries.

Conclusions
In the current study, we have presented experimental and theoretical results obtained for new Schiff base-type compounds. The reaction pathways and synthetic procedures have been presented. The physico-chemical properties of the compounds were characterized by NMR and IR spectroscopy. The biological activity was investigated based on diverse biological assays. It was found that Schiff base 13 inhibited the activity of both isoenzymes, COX-1 and COX-2 at a lower concentration than standard drugs, and its COX-2/COX-1 selectivity ratio was detected to be similar to meloxicam. The results of cytotoxicity assay showed that all of the tested compounds exhibited potent anti-cancer activity against A549, MCF-7, LoVo, and LoVo/Dx cell lines.
Finally, the quantum-chemical based DFT method simulations were performed. Two main possible conformations were investigated, and it was found that lower relative energy corresponded to the sandwich-like structure. As the last step of the study, the flexible docking method was applied and it was found that in comparison to known inhibitors of cyclooxygenases and affinity energies of docking meloxicam, piroxicam and rofecoxib to the receptors were similar to most of the energies obtained for the investigated structures.
These results indicate that the examined ligands could exhibit a significant activity towards human COX-1 and COX-2 receptors.
Considering the results obtained, we wish to underline that the combination of 1,3,4oxadiazole and hydrazide-hydrazone pharmacophore moieties has yielded promising results, and additional studies on these compounds are necessary. The research into their effects on apoptosis, the cell cycle and inflammation is planned to improve our understanding of their mechanisms of action.